Insulator-Superfluid transition of spin-1 bosons in an optical lattice in magnetic field 
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We study the insulator-superfluid transition of spin-1 bosons in an optical lattice in a uniform 
magnetic field. Based on a mean-field approximation we obtained a zero-temperature phase diagram. 
We found that depending on the particle number the transition for bosons with antiferromagnetic 
interaction may occur into different superfluid phases with spins aligned along or opposite to the 
field direction. This is qualitatively different from the field-free transition for which the mean-field 
theory predicts a unique (polar) superfluid state for any particle number. 



Q . I. INTRODUCTION 

The insulator-superfluid transition is an example of quantum phase transitions which take place at zero temperature. 
At the transition point the quantum ground state of the system changes in some fundamental way. This is accomplished 
by changing some parameter in the Hamiltonian of the system. Possible examples are transition between different 
quantized plateaus in a quantum-Hall effect under the change of magnetic field, transition between different phases in 
high-T c superconductor compounds under the variation of doping, metal-insulator transition in a conductor under the 
change of disorder, or phase transitions in baryonic matter under the change of density (e.g., transition from nuclear 
matter to a uniform liquid of neutrons, protons and electrons at the interface between a crust and a core of a neutron 
star). 

Originally the insulator-conductor transition has been studied in Fermi systems. However, analogous phenomena 
occur for a system of bosons, for example, in 4 He absorbed in porous media or granular superconductors in which 
the Cooper pairs may be considered as bosons. For pure Bose systems at zero temperature the conducting phase is 
the superfluid, so that the insulator-conductor transition corresponds to the onset of superfluidity. In the context of 
insulator-superfluid transition in liquid 4 He the problem has been studied in [1,2]. Recently such a problem attracted 
. much attention due to experimental realization of Bose-Einstein condensation (BEC) of atoms in magnetic and optical 
traps. It has been proposed that the insulator-superfluid transition might be observable when an ultracold gas of 
atoms with repulsive interactions is trapped in a periodic optical lattice [3]. Recently the transition in 87 Rb condensate 
• in a three-dimensional simple-cubic potential was realized experimentally by changing the potential depth [4]. An 
artificial optical lattice has a great advantage in that the system is basically defect-free which provides ideal conditions 
■ for study the quantum phase transitions. 

In the insulator phase the atoms are localized (by the effect of the potential) in the lattice sites so that fluctuations 
in the atom number on each site are suppressed. In such a state there is no phase coherence across the lattice. 
In addition the Mott insulator phase is characterized by the existence of an energy gap for the creation of particle 
or hole excitations, i.e., for the addition of particles to, or removal of particles from, the system. In a periodic 
potential the bosons can move from one lattice site to the next by quantum tunneling. The heights of the barriers 
can readily be adjusted since they are proportional to the intensity of the laser beam. The competition between 
the tunneling amplitude (kinetic energy) and the atom-atom interactions controls the point of the phase transition. 
When the tunneling coupling becomes large compared to the atom-atom interactions a delocalized wavefunction 
minimizes the dominant kinetic energy and the system undergoes a phase transition into the superfluid state. In 
such a state the atom number per site is not fixed and the wavefunction exhibits long-range phase coherence. Recent 
theoretical studies include a delocalizing transition of BEC in optical lattices [5] , the superfluid density and the quasi- 
momcntum distribution [6,7], the superfluid to Mott insulator transition with high filling factors (multiband situation) 
[8], commensurate-incommensurate transition [9], as well as Monte Carlo simulations for the Bose-Hubbard systems 
[10,11]. 

An advantage of an optical trap is that it liberates the spin degrees of freedom and makes possible condensation 
of spinor bosons. For example, 23 Na, 39 K and 87 Rb atoms (with nuclear spin / = 3/2 and electrons at s orbits) at 
low temperatures behave as simple bosons with a hyperfine spin / = 1. Optically trapped BEC creates a possibility 
to study a variety of phenomena in spinor many boson systems which posses extremely rich physics. The insulator- 
superfluid transition in spinor BEC is one of the possible novel phenomena. Recently the insulator-superfluid transition 
of spin-1 bosons interacting antiferromagnetically in an optical lattice was investigated theoretically [12,13], as well 
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as the transition in a multicomponcnt BEC system [14] and spin-2 Bose atoms [15]. In this paper we study the phase 
transition of spin— 1 bosons in the presence of an external magnetic field. Magnetic field lifts a degeneracy of the 
antiferromagnetic ground state and enriches the phase diagram. We show that magnetic field results in a qualitatively 
new effect. Depending on the field value the transition for bosons with antiferromagnetic intcrparticle interaction 
may occur into different superfluid phases with spins aligned along or opposite to the field direction. 



II. BASIC FORMALISM 

We consider a dilute gas of trapped bosonic atoms with hyperfine spin / = 1. The system is described by the 
following Hamiltonian [16]: H = H$ + Ha, where the spin symmetric and asymmetric parts are 

h s = Y1 I M t (-£f + v *) *« + yE / *tty*«*e<Pr, (l) 

a J ^ ' a,l3 J 

Ha = y J d * r *i~ + *±i*±i*-i*-i + 2*+*+*i* + 2*t 1 *+*_i* - 2*+*t 1 *i*_i + 

+2*+*+* 1 *_ 1 +2*+*t 1 * *o) , (2) 

here * Q (a = —1,0,1) is the atomic field annihilation operator associated with atoms in the hyperfine spin state 
|/ = l,m/ = a >, Vtr = Vo(sin 2 kx + sin 2 ky + sin 2 kz) is the optical lattice potential which is assumed to be the 
same for all three spin components, k — 2n/X, A is the wavelength of the laser light, Vo is the tunable depth of the 
potential well. The coefficients A s , A a are related to scattering lengths ao and a 2 of two colliding bosons with total 
angular momenta and 2 by A s = Airh 2 (a + 2a2)/3M and A a = 4-irh 2 (a2 — ao)/3M, where M is the mass of the 
atom [17]. For scattering of 23 Na ai = 54.7as and ao = 49.4as, where clb is the Bohr radius [18]. This suggests 
A a > (antiferromagnetic interaction). For 87 Rb a 2 = (107±4)as and ao = (110±4)as [17], that is A a < and the 
interaction is ferromagnetic. 

For single atoms the energy eigenstates are Bloch wave functions. An appropriate superposition of Bloch states 
yields a set of Wannier functions which are localized on the individual lattice sites. Expanding the field operators in 
the Wannier basis 

*«( r ) = ^a ai w(r - n), (3) 

i 

where a a i correspond to the bosonic annihilation operators on the ith lattice site, and keeping terms in first order in 
the hopping matrix element, the Hamiltonian (1), (2) reduces to the Bose-Hubbard Hamiltonian 

A s = - J E E + \sY. Mfti - 1). (4) 

<i,j> a i 



i 

(5) 

where Ni — a^i^ai is the number of atoms at lattice site i, < i,j > stands for the nearest-neighbor sites, 
2A Sj0 = A Si a J d 3 r\w\ 4 corresponds to the strength of the on-site interaction between atoms and 

J = - j dW(r-r t ) (-^+V t ^w(v-v j )>0 

is the hopping matrix element between adjacent sites This describes tunneling of atoms and is related to the 
kinetic energy. The parameter J exponentially depends on the strength of periodic potential Vo and can be varied 
experimentally by several orders of magnitude. One should mention that when number of particles in each well N is 
large enough the parameters A SiQ can be reduced by a factor iV~ 3 ' 5 . This occurs when condensates in each well reach 
the Thomas-Fermi regime. However, in experiments on 3D optical lattices N < 100 and the Thomas-Fermi regime is 
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usually not achieved. Higher filling factor can be realized in 2D and ID lattices [19,20]. Also, equation (3) suggests 
that atoms in different spin states are approximately described by the same coordinate wave function which is the 
case when the spin symmetric interaction is strong compared with the asymmetric part: |A S | 3> |A a | [16]. This is 
relevant to experimental conditions for 23 Na and 87 Rb atoms [18,17]. Introducing the operators 

L + = -\/2(ai~ao + dja_i), L_ = v^ajai + a^ao), L z = a^di — a^ft-i, L 2 = L + L_ + L 2 — L z 
which obey the angular momentum commutation relations [L + ,L_] = 2L Z , [L Z ,L±] = ±L± one can obtain [16] 

H A = \ a J2&-m). (6) 

i 

The operators Lj and L Z i commute with iV, and each other. 

To study the phase transitions, it is more convenient to perform calculations in the grand-canonical ensemble, so 
we add the term with the chemical potential fi. We also include an external uniform magnetic field B along the z 
axis [21]. Then the Hamiltonian becomes 

H = -J E + E [^M^ - 1) + A a (L 2 - 2N t ) - jj,Ni - bL zl ] , (7) 

<ij> a i 

where b = gfJ-BB which we assume to be positive, g is a Lande factor of an atom and \ib is the Bohr magneton. As 
we will see later the most interesting effects occur when the magnetic energy g/isB is less than or comparable with 
the energy of spin asymmetric interaction. For such values of magnetic field the energy of spin symmetric interaction 
is much larger than gusB and, hence, the description of atoms in different spin states by the same coordinate wave 
function is justified for such magnetic field. 

The consistent mean-field theory we will use corresponds to the following decomposition of the hopping terms 
[22-24]: 

& ai^j =< a+ L > a aj + a+ < a aj > - < a+ >< a aj >= ^ Q (a+ + a aj ) - ip 2 a , (8) 

where ip a =< >=< a a i > is the three component supcrfluid order parameter. In such a decomposition we omit 
the higher-order fluctuations (a^ — ^ a ){a a j — ip a )- Also we assume the order parameter to be real and neglect the 
Josephson-type tunneling term. Then Eq. (7) yields 



H = zJN s J2€+Y, 



\ s Ni(fii - 1) + A a(L 2 - 2Ni) - fj,Ni - bL zl - zjJ2Ma+ t + a a 



(9) 



where z = 2d is the number of nearest-neighbor sites, d is the space dimension and N 8 is the total number of lattice 
sites. This Hamiltonian is diagonal with respect to the site index i, so one can use an effective onsite Hamiltonian 
Hf 

Hf = - 1) + A a (L 2 - 2Ni) - iiNi - bL zi + z jJ2^ 2 a - zJ^M^ + (10) 

a a 

One should note that the simple mean-field approximation we use may not catch some interesting features discussed 
in [12] for a field-free transition. However, our simple approximation is good enough to describe a new effect we 
predict for a transition in magnetic field and, hence, there is no need to consider more complicated theories. 



III. PERTURBATION EXPANSION OVER THE SUPERFLUID ORDER PARAMETER 

We will treat the last two terms in Eq. (10) as a perturbation. For the second order phase transition the order 
parameter ip a continuously changes from zero (in insulator phase) to a finite value (in supcrfluid phase). Hence, in 
the vicinity of phase transition the order parameter is infinitesimally small and, therefore, to find the transition point 
the last two terms in Eq. (10) can be treated as a perturbation. In zero order approximation (we drop the site index 
i since the effective Hamiltonian is diagonal with respect to the lattice sites) 

Hf = X S N(N - 1) + A a (L 2 - 27V) - »N - bL z . (11) 
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The eigenstates of Hq S are states with defined total number of atoms in the site N (an integer number), the total spin 
per site I and its projection m along the magnetic field. We denote the eigenstates as \N, I, m >. The corresponding 
eigenvalues are 

E = X S N(N - 1) + X a [l(l + 1) - 2N] -nN- bm. (12) 

For b > the state with m = I corresponds to the lowest energy. We assume that all the atoms in the well are in 
the same orbital state, which is the ground state for the confining potential. Then the symmetry of the bosonic wave 
function and the structure of the operator L enforces a constraint on I. The allowable values of I are / = 0,2,4, ...N 
if N is even and I = 1, 3, 5, ...N if N is odd (also m = 0, ±1, ±2, . . . ,±l for any allowable I). This result is known in 
cavity QED; the details of the proof are provided in Appendix A of Ref. [25]. The eigcnfunctions \N,l,m > in terms 
of the Fock states \n\,no,ri-\ > with defined number of particles n a with the spin projection a are given by 

\N,l,m >= A k\k,N- 2k- m,k + m >, (13) 

k 

where A k satisfy the recursion relation 

A k+1 + l)(k + m+l){N -2k-m- l)(N -2k- m) + A k - iy /k{k + m)(N - 2k - m + 1)(N - 2k - m + 2) + 



+A k [k(N — 2fc — m + 1) + (N — 2k — m)(k + m+l)] = [1(1 + 1) - m(m + l)]A k /2. 
In some particular cases this equation has simple solutions. E.g., if m = I the coefficients are 




while for in = I — 1 



The coefficient Ao can be obtained from the normalization condition ^ fe A\ = 1. 

A. Mott ground state of spin-1 bosons 

The ground state of Hq S depends on the relation between A s , A a , b and /i. It is determined by minimizing the 
energy (12) with the constraint N + 1 = even. We assume A s > 0. When A a < the ground state is \N, N, N >, it is 
a ferromagnetic state in which all bosons occupy the m = 1 state. The number of particles N is determined from the 
condition 

2(N - 1)(A S + AJ < M + b < 2N(\ S + AJ. (14) 

When A a > the nature of the ground state depends on magnetic field. For 6 = the ground state is antifer- 
romagnetic. However, the number of particles per site depends on the relation between A s and A a . There are two 
possibilities: 

1) \ s > 2A a , then the particle number is determined from 

2(N - 1)A S - 4A Q < n < 2N\ S , N = 0,2,4,..., 1 = ( 
2(N - 1)A S < n < 2NX S - 4A a , N = 1,3,5,..., 1 = 1' { ' 

The number of particles at each site increases by +1 any time the chemical potential \i crosses the points 0, 2A S — 4A a , 
4A S , 6A S -4A a , ... 

2) A s < 2A a , then states with even N are the only possibility and 

(2N - 3)A S — 2A a < fi < (2N + 1)A S - 2A a , N = 0,2,4,..., 1 = 0. (16) 
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The particle number increases by +2 each time the chemical potential fi crosses the points A s — 2A a , 5A S — 2A a , 
9A;, — 2A a , 13 A s — 2A a , .... For A a > and 6 = the insulator-superfluid transition has been studied in detail in [13]. 
For A a , b > and given TV the number I is determined by the parity of N and the condition 

(2J-1)A„ <b< (2l + 3)X a . 

When b > A a (2iV — 1) the ground state is ferromagnetic with I = N. Let us define Z max so that A a (2Z max — 1) < b < 
A a (2Z max + 1). Then for b > the particle number is determined from (we assume A s > A a > 0) 

N < J max - 1 : 2{N - 1)(A S + A ) < /i + b < 2N(X S + X a ), 
N > ; max : 2(N - 1)A S - 2A a + (-1)"+'— (2Z max A a - b) < y. < 2NX S - 2X a - (-1) W+/ — (2Z max A a - b). 



B. Phase diagram 

We calculate the phase boundary between the insulator and the superfluid states for b > using second-order 
perturbation theory The idea is to find the energy as a function of the superfluid order parameter ip a : E = 
Eq + B\ipi + BqiJjq + B^iip 2 _ 1 + . . ., where coefficients B a depend on the system's parameters. In the insulator phase 
all B a > and, hence, ip a = minimizes the energy. Phase transition into a superfluid state occurs when one of 
the coefficients B a becomes negative. In Appendix A we calculate the necessary matrix elements of the perturbation 
operator between the ground state \N,l,l > and excited states of . Here we list the nonzero matrix elements: 



<N + 1,1 + 1, l\a£\N, I, I >= J N ^ l ^ 3 , < N — 1,1 + 1, l\a \N, I, I > 



21 + 3 



21 



<N+l,l + l,l + l\at\N,l,l>= 



{l + l)(N + l + 3) 



21 + 3 



< N- 1,1 + 1,1- l\ax\N,l,l >- 



N 



(2/ + 3)(2/ + l) : 



< N — 1,1 — 1,1 — lja.il JV, I, I >- 



l(N + l + l) 
21 + 1 



, <N + 1,1 + 1,1-11^^,1,1 >-- 



1 N + l + 3 
(2Z + 3)(2Z + 1) : 



<N + l,l- l,l-l\at 1 \N,l,l>= - 



l(N-l + 2) 
21 + 1 : 



< TV — 1, Z + 1, Z + l|a_i|7V, Z, Z >= 



(N-l)(l + l) 



21 + 3 



The second order perturbation theory results in the following expression for the energy correction to the \N,l,l > 
state in terms of the superfluid order parameter tp a 

E = X S N(N - 1) + X a [l(l + 1) - 2N] -nN-bl + zJ^ife- 



(2Z + 3)(2Z + 1) 



(I + 1)(N + I + 3)(2l + 1) N-l 

2NX S + 2lX a -b-fi + 2(1- N)X s + 2{l + 2)X a + b + fi + 2(1 - N)X S - 2(1 - l)X a + b + Ji 



l(N + 1 + 1)(2Z + 3) 



21 + 3 



N + l + 3 



N-l 



2NX s + 2lX a -fi 2(l-N)X s + 2(l + 2)X a +fi 



(zJ)H\ 



(2Z + 3)(2/ + l) 



^ + ^ + 3 



+ 



Z(A^ - Z + 2)(2Z + 3) 



(N-l)(l + l)(2l + l) 



2NX S + 2lX a + b - [i 2NX S — 2(1 + l)X a + b — fi 2(1 - N)X S + 2(1 + 2)X a - b + fj, 



(17) 



The leading correction to the energy is quadratic in ip a which demonstrates close analogy of the approach with 
Schrieffer- Wolff transformation. The idea of such canonical transformation is to eliminate the perturbation term in 
the Hamiltonian to first order [26]. Phase transition into a superfluid state with ip a ^ occurs when the coefficient 
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in front of ty\ becomes negative. For b = and A a > the superfluid transition occurs into a polar (spin— 0) state; 
the result is valid for both even and odd number of particles per site [13]. In such a transition the order parameter 
V>o becomes nonzero. 

When b > the situation changes substantially. One can see from Eq. (17) that as soon as the ground state 
corresponds to I ^ the coefficient in front of in the expression for E becomes smaller than the coefficient in front 
tpQ. At the same time the terms containing ipi and ip-i compete with each other. Depending on the value of b and 
the chemical potential (particle number) the system undergoes a transition into a superfluid state with Vi 7^ or 

The coefficient in front of tp 2 becomes zero when 

(2Z + 3)(2l + l) _ (Z + 1)(JV + Z + 3)(2Z + 1) N-l 

zJ ~ 2NX S + 2lX a - b - n + 2(1- N)X s + 2(l + 2)X a + b + fi + 

l(N + l + l)(2l + 3) 
2(l-N)X s -2(l-l)X a + b + ^ ( ' 

while for the coefficient in front of ip 2 _ 1 the condition is 

(21 + 3)(2Z + 1) (N + I + 3) (N — + 1)(2Z + 1) 



zJ 2NX S + 2lX a + b - fj, 2(l-N)X s + 2{l + 2)X a -b + fi 

l(N-l + 2)(2l + 3) 
+ 2NX s -2(l + l)X a + b- n' 



+ 



(19) 



The lowest value of J from Eqs. (18), (19) determines the point of phase transition. In Fig. 1 we plot the phase 
boundary in the J — /i plane for the insulator-supcrfluid transition in a magnetic field. In our estimates we take 
A Q = 0.4A S and b — 4.25A a . For J = (no tunneling) the system is in the ground state for a single well. The 
number of particles per site N as well as the total spin I depend on the chemical potential n. Solid line corresponds 
to a transition into superfluid phase with ipi ^ 0, while along the dash line the transition occurs into a state with 
■0-1 7^ 0. The phase boundaries between different superfluid phases with either non-zero ipi or are beyond the 
scope of the present paper and are not shown. In Fig. 2 we plot the phase diagram for A a = 0.4A S and b = 6A a . 
The possibility of the phase transition into different superfluid states can be understood as follows. In the insulator 
phase particle permutation symmetry imposes the constraint N + 1 = even; the parity of I is fixed by the parity of N. 
However, the interaction energy X a l(l + 1) — bl can be lower for I with the opposite parity. The symmetry constraint 
is relaxed in the superfluid phase. Due to the appearance of a coherent superfluid component the particle number 
per site is no longer fixed. Effectively the superfluidity removes the restriction on the Z's parity and, as a result, 
the superfluid order parameter may correspond to different projections of spins along the magnetic field in order to 
decrease the excessive interaction energy. The effect takes place for small magnetic field or large particle number per 
site when X a (2N — 1) > b. The coefficient A a can be estimated as A a = (7r 2 /3)i?_R[(a2 — a ) / A](Vb / 'Er) 3 / 4 , where 
En = 2ir 2 h 2 /MX 2 is the recoil energy [12]. For 23 Na atoms a 2 — ao ~ 5.3a_e and typical experimental parameters 
[27] A = 0.985^m, E R = 4.4 x 1CT 7 K, V ~ 12E R , N ss 3 we obtain the following limitation on the magnetic field 
B < 1CT 4 G. Since the current capability of magnetic shielding can reach 10~ 5 G, the low magnetic field regime is 
attainable. Experimentally the transition into different superfluid phases may be detected, for example, by creating 
a weak gradient of magnetic field. In such a case the superfluids with different magnetization will move in opposite 
directions. 

For b > X a (2N — 1) or A a < the ground state is ferromagnetic with I — N. In this case the coefficient in front of 
tp\ in the expression for the E is smallest among the coefficients in front of ip 2 , ip 2 , ip 2 _ 1 - This means that the order 
parameter ipi always appears first and the system undergoes a transition into ferromagnetic superfluid with ipi ^ 0. 
The transition occurs when 

,J W± + A a ) - b - ,i)(b + M - 2{N 1)(A S + A a )) Wfc+„^9/VQ r<?m 

z,J = — — , 2(7V - 1)(A S + A a ) < b + [i < 2N(X S + A a ). (20) 

Z(A S + A a J + + /i 

In Figs. 1 and 2 the ferromagnetic state is realized for small particle numbers per site (N < 2 and N < 3). 

One should mention that we neglected possible presence of slowly varying (at the length scale of the lattice period) 
component of the trapping potential. Such potential effectively provides a scan over fj, of the phase diagram at fixed 
values of J/X s and A a /A s . Also at higher value of the tunneling amplitude J additional phase transitions between 
different superfluid states may occur. Their description requires estimate of the next order terms in the energy and 
lies beyond the scope of our paper. One of the interesting question for future studies is superfluidity of "holes" which 
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might occur when some sites remain empty after filling the optical lattice by atoms. Also, thermal fluctuations may 
become important at nonzero temperature and destroy long-range phase coherence of superfluid. As a result, the 
transition can first occur into a conductor phase (with no long-range coherence) and only later into a superfluid state. 
This work was supported by NASA, Grant No. NAG8-1427. 

APPENDIX A: CALCULATION OF MATRIX ELEMENTS 

In the presence of magnetic field the ground state of zero-order Hamiltonian (11) is \N, I, I >. Let us consider matrix 
elements of the perturbation operator a+ + a Q between the ground state and the other eigenstates. To calculate matrix 
elements we consider the eigenstates of the Hamiltonian (11) in a form [28] 

\N, I, I >= ^ /( ^ = (fij-)' (e + ) {N - l)/2 \vac >, (Al) 
where + = &q 2 — 2a^a't 1 and / is a normalization factor 

States with lower magnetic quantum numbers can be obtained by operating L_ = ^/2{a^ax + a^ao), L + = 
V2(afa + a^a-i) to \N,l,l > and using L-\N,l,m >= ^J{1 + m){l - to + Tj\N, I, m - 1 >, L + \N,l,m >= 
\J{1 + m + — to) | A, I, m + 1 >. Another useful operator is L z = a^fii — aj^a-i, L z \N,l,m >= m\N 7 l 7 m >. 
We note the properties 

[L±,e+]=0, [4,6+]=0, [L_,(a+) 1 ] = V2la+(a+) l -\ [L_,a+] = ^at 1 , 



6+|A, l, l > = y/(N - I + 2){N + l + 3)\N + 2,l,l>, L±\vac>= 0, L z \vac>= 0, 

[oi, (6+)"] = -2nal 1 (6+)"- 1 , [a , (0+)"] = 2na+(6+)"- 1 , [o_i, (6+)"] - -2na+(e+)"-\ 

[L + ,a± 1 ] = V2a£, [L+, ai] = -V2a , [L 2 , alj = -a^, ai] = -o_i, 

where [a, 6] stands for the commutator of a and b. Then using (Al) we obtain 

y/f(N-l,l-l)a+\N - 1, 1 - 1, Z - 1 >= a+(a+)'- 1 (e + ) (A '-' )/2 |mc >= ^=-[£-> (ai")'](& + ) (JV_,)/2 |vac >= 



^L_(a+) i (e+)( N - ( )/ 2 |m C >^ y^Rt_\N,l,l >= ^ fi * J) \N,l,l - 1 > 



or 



Further 



a+\N,l,l>=J N +^ + 3 3 \N+l,l + l,l>. (A2) 



a+|AT, M >= 1 (q|)'^(9+)(^')/> flC >= J m +l'l +1) 

v r iWJ) V w) 



JV+l,i + l,f + l >= 



" + 1 ^V + 3 W U + U + 1>. (A3) 
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a_,\N, I, I >= 1 (o+)'a 1 (A+)( jV -')/ 2 |„». >= - ( A " l 2 (»+)'+ 1 (A+)( jV -'- 2 )/ 2 |„», >= 



a |7V, i, J >= 1 (g+)'a (e+)( jv -')/ 2 | mc >= ( A "~ (g+)'g+(e+)( jv -'- 2 )/ 2 | mc > = 



= {N-^i^^^N-2,l,l>=^N-l,l + l,l>. (A5) 
Now let us find a't 1 \N, 1,1 >. Using 

V2a£\N,l,l >= L+at^NJJ >, 

we obtain 



a± 1 \N,l,l>=^ j^^±^\N + 1,1 + 1,1-1 >+B\N +1,1 -1,1-1 >, (A6) 

where B is a coefficient. Other terms such as \N + 1, m, m > do not enter Eq. (A6). This can be proven by applying 
L z to Eq. (A6) and using the property L z a't 1 \N,l,l >= (I — l)al 1 |A r , 1,1 >. To find the coefficient B one can 
apply a_i to Eq. (A6) and consider the matrix element < N, I, Z|a_ia! 1 |./V, I, I >= 1+ < N, I, Z|a^ 1 a_i|A'', I, I >= 
1 + (N - + I)/ {21 + 3). As a result, we obtain B = -y/l(N - I + 2)/(2Z + 1). Hence 



at.l^, / >= ^^+^^1^ + 1, 1 + 1, / - 1 > - ^^(^+2) , ^ + 1, Z - 1, 1 - 1 > . (A7) 
Finally let us calculate d\\N,l,l >. Using 

-V2a \N,l,l >= L + ai\N,l,l >, 

we find 



/ N - I 

ai\N,l,l>=-J \N -\,l + \,l-\> +C\N - \,l-\,l-\> . 

The coefficient C can be obtained by applying af and considering the matrix element < N,l,l\didi\N,l,l >. Ulti- 
mately 



&1 '^ >= -V (2Z + ^ (A8) 
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FIG. 2. The same as in Fig. 1, but for b = 6A a and A a = 0.4A, 
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